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Abstract 

The connection between analytic and Monte Carlo calculations of soft gluon 
emission is reanalyzed in light of recent, theoretical developments in resum- 
'nI" ' mation. An alternative Monte Carlo algorithm is suggested which incorpo- 

^^ . rates (1) corrections beyond leading order to the showering and (2) smoothly 

ON ■ 

Q> ■ merges with the higher order calculation of single, hard parton emission. In 

i-C I particular, it is possible to study jet properties in heavy boson production 

P-i! 

for all values of the boson transverse momentum with a total cross section of 



NLO accuracy. The specific cases of W, Z and Higgs boson production at the 
Tevatron are addressed using a modified version of the PYTHIA Monte Carlo. 
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I. INTRODUCTION 

In the near future, experiments at two, high-energy hadron colhders will search for 
evidence of physics that supersedes the standard model. Important among the tools that 
will be used in these searches are showering event generators or showering Monte Carlos 
(SMC's). Among the most versatile and popular of these are the Monte Carlos HERWIG [|l|, 
ISA JET 0, and PYTHIA 0. SMC's are useful because they accurately describe the emission 
of multiple soft gluons, which is, in effect, an all orders problem in QCD. However, they 
only predict total cross sections to a leading order accuracy, and, thus, can demonstrate a 
sizeable dependence on the choice of scale used for the parton distribution functions (PDF's) 
or coupling constants (particularly a^). Also, in general, they do not translate smoothly into 
kinematic configurations where only one, hard parton is emitted. In distinction to SMC's are 
certain analytic calculations which account for multiple soft gluon emission and higher order 
corrections to the hard scattering. These resummation calculations, however, integrate out 
the kinematics of the soft gluons, and, thus, are limited in their predictive power. They can, 
for example, describe the kinematics of a heavy gauge boson produced in hadron collision, 
but cannot predict the number or distribution of jets that accompany it. However, searches 
for new physics, either directly or indirectly through careful measurements of standard model 
predictions, often demand detailed knowledge of kinematic distributions and jet activity. 
Furthermore, VT+jets (and Z+jets) processes are often backgrounds to SUSY or technicolor 
signatures, and we demand a reliable prediction of their properties. The aim of this present 
study is to show how the positive features of the analytic resummation calculations can be 
used to improve the showering algorithms. 

The outline of the remainder of this paper is as follows. First, we review the parton 
shower (Sec. II) and the analytic resummation formalisms (Sec. HI). We then show how 
to modify the showing algorithm to incorporate higher-order corrections to the total cross 
section (Sec. IV). Furthermore, we show how to correct the soft gluon approximation made 
in the showering, so that there is a smooth transition between the showering and hard 
emission limits (Sec. V), and we compare this approach to other work. Finally, we present 
numerical results (Sec. VI) for W, Z and Higgs boson production at the Tevatron in Run 
II using a modified version of PYTHIA, and our conclusions (Sec. VII). 

In the ensuing discussion, we focus on the specific example of W boson production at 
a hadron collider, when the W decays leptonically. The results apply equally well to '~f*, Z 
and Higgs bosons (or any heavy, color-singlet particle) produced in hadron collisions. 



II. PARTON SHOWERS 

SMC's are based on the factorization theorem [Q, which, roughly, states that physical 
observables in any sensible gauge theory are the product of short-distance functions and 
long-distance functions. The short-distance functions are calculable in perturbation theory. 
The long-distance functions are fit at a scale, but their evolution to any other scale is also 
calculable in perturbation theory. 

A standard application of the factorization theorem is to describe W boson production at 
a pp collider at a fixed order in as- The production cross section is obtained by convoluting 
the partonic subprocesses evaluated at the scale Q with the PDF's evaluated at Q. The 
partons involved in the hard collision must be sufficiently virtual to be resolved inside the 
proton, and a natural choice for the scale Q is Q = My/ 0. Prior to the hard collision, 
however, the partons are not resolvable in the proton (i.e. the proton is intact) and have 
virtualities at a much lower scale Qo of the order of 1 GeV. The connection between the par- 
tons at the low scale Qo and those at the high scale Q is described by the DGLAP evolution 
equations 0. The DGLAP equations include the most important kinematic configurations 
of the splittings a —>■ be, where a, b and c represent different types of partons in the hadron 
{q,g, etc.). Starting from a measurement of the PDF's at a low scale Qq, a solution of 
the DGLAP equations yields the PDF's at the hard scale Q. Equivalently, starting with a 
parton c involved in a hard collision, it is also possible to determine probabilistically which 
splittings generated c. In the process of de-evolving parton c back to the valence quarks 
in the proton, a number of spectator partons (e.g. parton b in the branching a -^ be) are 
resolved. These partons constitute a shower of soft jets that accompany the ly-boson, and 
influence its kinematics. 

The shower described above occurs with unit probability and does not change the total 
cross section for PF-boson production calculated at the scale Q ^. The showering can be 
attached to the hard-scattering process based on a probability distribution after the hard 
scattering has been selected. Once kinematic cuts are applied, the transverse momentum 
and rapidity of the W^-boson populate regions never accessed by the differential partonic 
cross section calculated at a flxed order. This is consistent, since the flxed-order calculation 
was inclusive (i.e., pp -^ W + X) and was never intended to describe the detailed kinematics 
of the ly-boson. The parton shower, in effect, resolves the structure of the inclusive state of 
partons denoted as X. In practice, the flxed order partonic cross section (without showering) 
can still be used to describe properties of the decay leptons as long as the measurable is well 
deflned (e.g., the number of leptons with central rapidity and high transverse momentum. 



but not the distribution of rapidity or transverse momentum of the W). 

Here, we review parton showering schematically. More details can be found in Ref. [Q. 
First, for simplicity, consider the case of final state or forward showering, where the parton 
virtuality Q evolves forward to the low scale Qo- The basis for developing a probabilistic 
picture of final state showering is the DGLAP equation for the fragmentation functions: 
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where Pa~^bc is an unregularized splitting function, aatc is the coupling times color factor, 
and e is a cutoff. The equation can be rewritten as 

^^ (D,(x,g)/A(g)) = £ ^ ''^'f^'^^ Pa^Uz){Dbix/z,Q)/AiQ)) 

or, after integrating both sides of the expression, 

D^{x,t') = D.(x,t)A(t') + f f df- ^12 ""'f' ^"^ Pa^U^)Db{x/z, t"), (2.2) 
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where t = InQ^, with similar definitions for t' and t". The function 

A(t') = exp (- £ J'~' dt"dz '''''f^^"h ,^U^)] (2.3) 

gives the probability of evolving from the scale Q' = e* to Qq = e^° with no resolvable 
branchings, and is called the Sudakov form factor, to is a cutoff scale for the showering. 
A(t') contains all the information necessary to reconstruct a shower, since it encodes the 
change in virtuality of a parton until a resolvable showering occurs. Showering is reduced to 
iterative solutions of the equation r = A(t')/A(t"), where r is a random number uniformly 
distributed in the interval [0, 1], until a solution for Q' is found which is below a cutoff. For 
consistency, the cutoff should represent the lowest scale of resolvable emission Qq. 

For the case of initial state radiation, several modifications are necessary. The fragmen- 
tation function is replaced by a parton distribution function, and the evolution proceeds 
backwards from a large, negative scale — IQ^I to a small, negative cutoff scale — IQqI- There 
are two equivalent formulations of backwards showering based on the probabilities 

exp (- f r^ dt''dz^^^^^^^P,^U^)^^^] , x' = x/z [9], (2.4) 
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and 
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After choosing the change in virtuahty, a particular backwards branching is selected 
from the probability function based on their relative weights (a summation over all possible 
branchings a ^ 6c is implied these expressions), and the splitting variable is chosen by 
solving the equation 

fx/x' Ay ^ f\-e Ay ^ 

/ ^P_,,(z)/(x/z, t') = r / -P„^,,(z)/(x/z, t'), (2.6) 
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where r is a random number. The details of how a full shower is reconstructed in the 
PYTHIA Monte Carlo, for example, can be found in Ref. |^. The structure of the shower can 
be complex: the transverse momentum of the W^-boson is built up from the whole series of 
splittings and boosts, and is known only at the end of the shower, after the final boost. 

The SMC formulation outlined above is fairly independent of the hard scattering process 
considered. Only the initial choice of partons and possibly the high scale differs. Therefore, 
this formalism can be applied universally to many different scattering problems. In effect, 
soft gluons are not sensitive to the specifics of the hard scattering, only the color charge of 
the incoming partons. 

III. ANALYTIC RESUMMATION 

At hadron colliders, the partonic cross sections can receive substantial corrections at 
higher orders in a^. This affects not only the total production rate, but also the kinematics 
of the W boson. At leading order (a°), the VT-boson has a ^(Q^) distribution in Q'^. 
At next-to-leading order, the real emission of a single gluon generates a contribution to 
da /dQ\ that behaves as Qx^as{Q\) and Q^'^asiQ'r) ^"^{Q"^ IQt) while the leading order, soft, 
and virtual corrections are proportional to —5{Q\). At higher orders, the most singular 
terms follow the pattern of as(<5r)"E^"=o ln'"(gVg|.) = a^L™ = V". The logarithms 
arise from the incomplete cancellation of the virtual and real QCD corrections, but this 
cancellation becomes complete for the integrated spectrum, where the real gluon can become 
arbitrarily soft and/or coUinear to other partons. The pattern of singular terms suggest that 
perturbation theory should be performed in powers of V^ instead of a". This reorganization 
of the perturbative series is called resummation. 



The first studies of soft gluon emission resummed the leading logarithms |11,12|, leading 



to a suppression of the cross section at small Qt- The suppression underlies the importance 



of including sub-leading logarithms [|I^]. The most rigorous approach to the problem of 



multiple gluon emission is the CoUins-Soper-Sterman (CSS) formalism for transverse mo- 
mentum resummation [|T^, which resums all of the important logarithms. This is achieved 
after a Fourier transformation with respect to Qt in the variable 6, so that the series involv- 
ing the delta function and terms V^ simplifies to the form of an exponential. Hence, the soft 
gluon emission is resummed or exponentiated in this 6-space formalism. To be more correct, 
the Fourier transformation is the result of expressing the transverse-momentum conserving 
delta functions 5^'^\Qt — Z^^tJ in their Fourier representation. Also, the exponentiation 
is accomplished through the application of the renormalization group equation, not by re- 
organizing an infinite sum of terms. Despite the successes of the 6-space formalism, there 
are several drawbacks. Most notable for the present study is that it integrates out the soft 
gluon dynamics and does not have a simple physical interpretation. 

The CSS formalism was used by its authors to predict both the total cross section to 



NLO and the kinematic distributions of the PF-boson to all orders [T^ at hadron colliders. 



A similar treatment was presented using the AEGM formalism |]T6|, that does not involve 
a Fourier transform, but is evaluated directly in transverse momentum Qt space. When 
evaluated at NLO, the two formalisms are equivalent to NNNL order in a^, and agree 
with the fixed order calculation of the total cross section |T^. A more detailed numerical 



comparison of the two predictions can be found in Ref. |T8 



Recently, the AEGM formalism has been reinvestigated, and an approximation to the 
6-space formalism has been developed in Qr^space which retains its predictive features |T^ 
(see also the recent eprint [^). This formulation o?oe5 have a simple, physical interpretation, 
and can be used to develop an alternate algorithm for parton showering which includes 
higher-order corrections to the hard scattering. 

In the 6-space formalism, the differential cross section of the ly-boson produced in 
association with soft gluons is: 

da{hih2 -^ 1/(*)X) 1 



d'hd'-'^^Wih, Q, xi, X2) + Y{Qt. Q, xi, X2). 



(3.1) 
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In this expression, Q, Qt and y describe the kinematics of the boson V , the function Y is 
regular as Qt -^ and corrects for the soft gluon approximation, and the function W has 
the form: 
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and 



fe ® fl/h^ (xi,^) = /' ^QK^,C7i,C2,/i = Cs/b)fi/,,{^,,f, = C,/b). (3.4) 

In these expressions, Ci, C2 and C3 are constants, iif is a function that describes the hard 
scattering, and A, B, and C are calculated perturbatively in powers of as'- 

{A,B,C) = f:(^^y{A,B,Cr^ 



n=0 \ ^ 



(the first non-zero terms in the expansion of A and B are for n = 1). The functions C*^"^-* 
are the Wilson coefficients, and are responsible for the change in the total production cross 
section at higher orders. In fact, (C ® /) is simply a redefinition of the parton distribution 
function obtained by convoluting the standard ones with an ultraviolet-safe function. 

We remove the constants Ci, C2 and C3 from these expressions by choosing their canonical 
values |]n[, which also removes large logarithms. At leading order, the expression for the 



production of an on-shell Pi^-boson simplifies considerably to: 

da{hih2 -^ WX) ^ r (fb iiQ^ „g(^Q) f{xi,b)f{x2,b) 

dQi ^V (2vr)2' /(xi,g)/(x2,g)' ^^-'^ 

aQ = K j — -f{xi,Q)f{x2,Q), 
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where k contains physical constants and we ignore the function Y for now. The expression 
contains two factors, the total cross section at leading order o"o, and a cumulative probability 
function in Q"^ that describes the transverse momentum of the PF-boson (the total integral 
over Q"^ transforms e*'''*^^ to 5^'^\b)). Except for the complication of the Fourier transform, 
the term e-'^/^f{x,b)/f{x,Q) is analogous to A{Q)f{x,Q')/A{Q')f{x,Q) of the SMC. 

Equation (3.1) in 6-space has a similar structure in Qr^space. This is surprising, since 
the 6-space result depends critically on the conservation of total transverse momentum. To 
NNNL accuracy, however, the Qt space expression agrees exactly with the 6-space expres- 
sion, and has the form |jl9[: 

"'".Q^.Qj.'^Jr' = ^'^W-'?-^-"^^' + mr,Q,....). (3.6) 

Again ignoring Y, we can rewrite this expression as: 
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The factor ai is the total cross section to a fixed order, while the rest of the function yields 
the probability that the ly-boson has a transverse momentum Q^. 



IV. A MODIFIED SHOWERING ALGORITHM 

The primary result of this paper is to exploit the expression for the differential cross 
section in Eq. (3.7), which has the form of a leading order cross section times a backwards 
evolution. We generalize the function A{t)/f{x,t) x f[x,t')/A{t') of the standard back- 
wards showering algorithm to yW (the square root appears because we are considering the 
evolution of each parton line individually). 

To implement this modification in a numerical program, like PYTHIA, we need to provide 
the new, modified PDF (mPDF) based on the Wilson coefficients. At leading order, the only 
Wilson coefficient is C^j = 6ij6{l — z), and we reproduce exactly the standard showering 
formulation. For VT-boson production at NLO, the Wilson coefficients C are: 

C^l^ = S^k [1(1 -z) + li^r' - 8)6(1 - z)} , C« = i^(l - z). (4.1) 

To NLO, the convolution integrals become: 

^ a,(/i)l^_2 ^A 
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77(1 - z)fi{x/z,fi) + -(1 - z)fg{x/z,fi) 



(4.2) 
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and fg{x,fi) is unchanged. The first term gives the contribution of an unevolved parton to 
the hard scattering, while the other two contain contributions from quarks and gluons with 
higher momentum fractions that split q —* qg and g —>■ qq, respectively. 

We are relying on the fact that the Sudakov form factor used in the analytic expressions 
and in the SMC are equivalent. In fact, the integration over the quark splitting function in 
A{Q) yields an expression similar to the analytic Sudakov: 



f '"'dzCpl--^)=CF[\n\ ^1 - 3/2(1 -2^J~A«ln(gVg|) + S«, (4.3) 
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where Zm = (qT^ \ is an infrared cutoff, terms of order Zm and higher are neglected, and 



{Q+Qt) 



the z dependence of the running coupling has been ignored ||2T| . Note that the coefficients 
A^^"* (Cf) and 5'-^^ {—3/2Cf) are universal to qq annihilation into a color singlet object, just 
as the showering Sudakov form factor only knows about the partons and not the details of 
the hard scattering. For gg fusion, only the coefficient A^^'^ (3) is universal. In general, at 
higher orders, the analytic Sudakov is sensitive to the exact hard scattering process. 

While the Sudakov form factors are similar, there is no one-to-one correspondence. 
First, the Q^^space Sudakov form factor is expressed directly in terms of the Qt of the 
heavy boson, while, in the SMC's, the final Qt is built up from a series of branchings. 
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Secondly, the integral on the left of Eq. (4.3) is positive (provided that Zm < \), while the 
analytic expression on the right can become negative. This is disturbing, since it means 
subleading logarithms (proportional to B) are dominating leading ones. In the exact SMC 
Sudakov, the kinematic constraints guarantee that A((5) < 1. In this sense, the Sudakov in 
the SMC is a more exact implementation of the analytic one. We feel that the agreement 
apparent in Eq. (4.5) is compelling enough to proceed assuming the two Sudakov form 
factors are equivalent. In our phenomenological analysis, we calculate A((5) numerically. 

V. HARD EMISSION CORRECTIONS 

The SMC and resummation formalisms are optimized to deal with kinematic configu- 
rations that have logarithmic enhancements L. For large Qt — Q, there are no such en- 
hancements, and a fixed order calculation yields the most accurate predictions. The region 
of medium Qt, however, is not suited to either particular expansion, in a^L"^ or a". 

The problem becomes acute for SMC's. In the standard implementation of SMC's, the 
highest Qt is set by the maximum virtuality allowed, Q = My/ in our example, so that 
the region Qt > Q is never accessed. However, at Qt > Q-, the fixed order calculation is 
preferred and yields a non-zero result, so there is a discontinuity between the two predictions. 
Clearly, the SMC underestimates the gluon radiation well before Qt — Q, but the fixed order 
calculation make equally overestimate the radiation if extended to the region Qt < Q- 

In the 6-space calculation, the same sort of behavior would occur, except that contribu- 
tions to the cross section that are not logarithmically enhanced as Qt -^ can be added 
back order-by-order in a.,. This procedure corrects for the approximations made in deriving 
the exponentiation of soft gluon emission. This correction is denoted Y, or Yf in contrast 
to the resummed term Yr. If the coefficients A and B are calculated to high-enough accu- 
racy, one sees a relatively smooth transition between Eq. (3.1) and the NLO prediction at 
Qt = Q- The situation is even better in the Qr^space calculation, since the matching at 
Qt = Q is guaranteed at any order. 

It is useful to review Y in the resummed calculation, which has the form 
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with 



A-..(^) = C^ (^) , P,^,{^) = l[z' + {l- zr] . (5.2) 

The first term in each expression is proportional to the squared matrix element for the hard 
emission, while the terms proportional to Q^^ are the asymptotic pieces from the Sudakov 
form factor. Similar corrections can be derived for the SMC. In general, the hard emission 
contributes to the hadronic cross section proportional to 

fa/hi{Xa,Q)fb/h2{^b,Q) O"0 ^exact, 

where (Tq is the leading order partonic cross section, and the Ag.xact are the expressions 
multiplying the delta functions in Eq. (5.2). The showering contribution to the same order 
(assuming showering of only one parton with virtuality —t) is 

dXadr]b6{xa'nb-Q^/S)fa/hi{xa,t)fb/h2{Vb,t) (^0 -^P{z)-^^p—^ — ^A(-t), 

tz fb/h2{Vb,q) 

or, by changing variables to Xb = Tjb/z, 

Z fa/hi{Xa,Q)fb/h2{Xb,Q) (^0 ^shower- 

As expressed in the functions R in Eq. (5.2), the single, hard emissions generated by the 
showering can be subtracted from the exact squared amplitude to include the remaining 
NLO corrections not present in the modified PDF. This defines a Y term for showering. 

The effects of Y are included by generating the NLO subprocesses and performing the 
subtraction of showering contributions on an event-by-event basis. We illustrate this ex- 
plicitly for the qq' — > Wg subprocess. Each event receives an additional weight fy before it 
is accepted or rejected (this is accomplished in the subroutine PYEVWT already provided in 
the PYTHIA code for the user to reweight any process), where fy is defined as: 

^exact ^shower 

\ ' 

^exact 



h - 4 
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Ae.act = ^^ ^'^'+j" ^'^' , A,,„^,, = P^^^{z) i—M-i) + ^^{-U)\ (5.3) 

tu \—t —u J 

with z = Q'^/s. The correction for qg -^ Wq' has a similar form: 

_ (g + t)^ + (t + n)^ _pr^^A^M rc;/!^ 

^exact — TT ; ^sfcoiuer — -Tg^g^Z) -/\[—U). l'3.4j 

One can show that the showering corrections are either smaller than the exact squared 
amplitude or equal to it in the limits when u or t —>■ 0. In writing this expression, we are 
ignoring the possibility that later emissions are harder than the first one [^. It has been 



shown previously that virtuality-ordered showering, as implemented in PYTHIA, yields the 



hardest emission first 90% of the time |27], and we ignore this technical detail. 

At this point, it is useful to compare the scheme outlined above to other approaches 
at improving the showering algorithm. One class of corrections is based on phase-space 
splitting, where part of a NLO matrix element is treated with LO kinematics and part with 



exclusive NLO kinematics p4| , p5| . The idea is to allow parton showering for one of these 
configurations, and not the other, but particular care must be taken not to mix the different 
regions of phase space. There is an adjustable parameter that splits the phase space, and 
physical observables are sensitive to the exact choice (see the discussion in Ref. ^^ regarding 



Qx^). In the approach of Ref. [24], the splitting parameter is tuned so that the contribution 
with LO kinematics vanishes. The resultant showering of the term with exclusive NLO 
kinematics generates emissions which are harder than the first "hard" emission, which is not 
consistent. Furthermore, the sphtting parameters must be retuned for different processes 
and different colhders. This scheme is guaranteed to give the NLO cross section before cuts, 
but does not necessarily generate the correct kinematics. 

The other class of corrections modifies the showering to reproduce the hard emission limit 
P^ , P7[] . While this can be accomplished, it does so at the expense of transferring events from 
low Qt to high Qt- There is no attempt to predict the absolute event rate, but only to 
generate the correct event shapes. In some implementations, the high scale of the showering 
is increased to the maximum virtuality allowed by the collider energy. This is contrary to 
the analytic calculations, where the scale Q = Mw, for example, appears naturally (in the 
choice of constants Ci, C2 and C3 which eliminate potentially large logarithms). This scheme 
will generate the correct hard limit, but will not generate the correct cross section in the 
soft limit. 

The present formulation contains the positive features of both schemes. Phase space is 
split (but without any adjustable parameters - which is also true for the corrections outlined 



in Ref. p3[), but higher-order corrections are also applied to the showering algorithm. This is 
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contained in the modified PDF. Instead of applying corrections to tlie sliowering to reproduce 
tlie liard limit, corrections are applied to the hard emission cross section to avoid double 
counting with the showering. The corrections allow a smooth transition between the explicit 
hard emission and the showering. 

So far, the discussion has been theoretical. In the next section, we demonstrate the 
phenomenological implementation of these ideas. 

VI. NUMERICAL RESULTS 

We have applied the showering modifications outlined above to W and Higgs boson pro- 
duction in Run I and Run II at the Tevatron, using a modified version of the PYTHIA Monte 
Carlo 0. Most of the previous discussion applies exactly to the case of Z boson produc- 
tion. In particular, the modified PDF used for the showering and the corrections /y are 
exactly the same. For Higgs boson production, we use the expressions for the Wilson coef- 
ficients presented in Ref. [^. Some technical points should be noted. First, we have not 



attempted to modify the Sudakov form factor implemented in PYTHIA; we have only modi- 
fied the PDF's that drive the showering. The generalization of showering to NLO has only 
been accomplished for final state showering |2^, and is technically complicated. Secondly, 
the expression for fy, with z = Q'^/s and A{Q) calculated numerically using Eq. (2.3), does 
not vanish fast enough as Qt —^ 0. This means that /y begins to increase sharply before 
vanishing at Qt — because of phase space. We do not believe this rise is physical, so we 
force the vanishing of /y below a certain cutoff. For all cases considered in this study, a 10 
GeV cutoff is adequate. We allow Y to shower with a maximum virtuality fixed at —t or 
—u, and include primordial transverse momentum for the incoming partons, so this cut is 
smeared out. Ideally, however, no cut would be necessary, and this issue is left for future 
study. Finally, we have only considered leptonic decays of the W and the 77 decay of the 
Higgs boson (to avoid any effects of final state radiation), but our final results are scaled to 
the total production cross section. Our results for Z boson production are similar to those 
for W boson production, so we do not comment on them further. 

For our numerical results, we present the Qx distribution of the heavy boson produced 
at the Tevatron. These distributions are in good agreement with analytic calculations, but 
cannot be predicted accurately by the standard showering algorithm. Secondly, we present 
jet properties for the same processes, which are not significantly altered from the predictions 
of the standard showering algorithm. These cannot be predicted by analytic calculations. 
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A. Heavy Boson Properties 

Our first goal is to test the predictions of our sliowering algoritlim on Run I data. In 
Fig. 1, tfie transverse momentum of tlie W boson (solid line) is shown in comparison to 
D0 data |]^ (the three dashed lines show the data and the upper and lower error esti- 
mates). The modified PDF (mPDF) is calculated using CTEQ4M PDF's. As in analytic 
calculations, the position of the peak from the SMC depends on non-perturbative physics 
1 3^. In PYTHIA, this is implemented through a Gaussian smearing of the transverse momen- 



tum of the incoming partons. To generate this plot, we have changed the default Gaussian 
width from .44 GeV to 3.6 GeV, which is more in accord with other analyses. This is the 
value used in all subsequent results. The numerical agreement between our prediction and 
the data is very good, except, perhaps, in the region of 30-50 GeV.Q Our prediction in this 
region depends on the details of the correction /y, and any excess is probably related to the 
bad behavior of fy at very small Qt- 

Next, we present our results on the Qt distributions for W and Higgs bosons {ttih = 100 
GeV) produced at the Tevatron in Run II. Figure 2 displays the Qt distribution of the 
W boson at Run II generated using the modified version of PYTHIA. The solid line is the 
full distribution, and the dashed and dotted lines show the individual contributions from 
the corrected showering and the corrected hard emission {Y) piece. For reference, we show 
the Qt distribution using the default version of PYTHIA and CTEQ4L PDF's. Because the 
maximum virtuality of the showering is set to the scale Mw, the SMC contributions are 
suppressed beyond Qt — Mw- Note how the Y piece fills in the intermediate Qt region 
down to small Qt, where showering gives the preferred result. The total cross section 
predicted by PYTHIA is (18.2, 20.5, 23.9) nb using CTEQ4L, CTEQ4M, and mPDF+F. The 
total increase in rate from LO to NLO is in good agreement with Ref. |p2| at v5 = 2 TeV. 



Our numbers for v5 = 1.8 TeV (21.2 nb using mPDF+y) also agree with the CDF and 
D0 data g. 

It is interesting that the mPDF calculation, without Y, yields nearly the same rate as 
using just the CTEQ4M PDF (21.6 nb vs. 20.5 nb). This is anticipated by the smallness 
of the virtual correction in the Wilson coefficients oc tt^ — 8. However, this need not be the 
case for different processes or different colliders. 



^This agreement also relies on using PYTHIA v6.125 or higher, which treats the showering kinemat- 
ics more correctly. This correction has an even greater impact on Higgs boson production, where, 
previously, the transverse momentum distribution was much broader than for W boson production. 
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The results for the production of a Higgs boson with mass rriH = 100 GeV is shown in 
Fig. 3. Here, the correction to the lowest order process is quite large. The total cross section 
predicted by PYTHIA is (0.50, 0.48, 0.94) pb using CTEQ4L, CTEQ4M, and the mPDF+F. 
The final number are in good agreement with Ref . [^ . We have used a primordial kr tuned 
to Run I data. However, if we believe that the non-perturbative function should have the 
same form as the perturbative Sudakov form factor, then the primordial kj- should scale like 
Ca/Cf relative to the qq case [^, and the peak of the Qt distribution would shift to the 
right. 

It is interesting to know if the kinematic distributions for the heavy bosons can be 
reproduced using the standard showering algorithms with a multiplicative ii'-factor that 
yields the total NLO rate. Figure 4 shows the ratio of the Higgs boson transverse momentum 
distributions calculated from mPDF+F to CTEQ4L times K {K > 2). There are variations 
as large as 10% in the important regions of small and medium Qt- Of course, the effect is 
much larger for the large Qt region where there is almost no rate from the standard parton 
showering. 

B. Jet properties 

In Figs. 5, 6 and 7, we present jet properties for the W and Higgs boson production 
processes. Jets are defined using the cone clustering algorithm of the PYTHIA subroutine 
PYCELL with a cone of size R = 0.5, Et > 5 GeV, and |?7| < 2.5. For W boson production in 
Fig. 5, we present the jet multiphcity distributions for all Qt, Qt < 10 GeV, 10 < Qt < 20 
GeV, and Qt < 30 GeV. We compare the showering predictions using mPDF (solid line), 
CTEQ4M (dashed line), and CTEQ4L (dotted line). To study the effects of the modified 
showering algorithm, we do not include Y, which would increase the Njet = 1 bin of the 
mPDF prediction. All distributions are normalized to unity. From Fig. 5, we see that the 
predictions have only minor differences, which is expected since the Wilson coefficients for 
W production are nearly unity. 

For the case of Higgs boson production in Fig. 6, we study the regions Qt < 25 GeV, 
25 < Qt < 50 GeV, and Qt < 75 GeV. There are more noticeable differences, and much 
more radiation in general than for the W boson case. The higher order PDF's generate 
slightly more jet activity and yield similar distributions. In general, though, there are no 
dramatic changes in the distributions. This is not too surprising, since the showering depends 
on the ratio of the modified PDF's evaluated at two different scales, which is not as sensitive 
to the overall normalization of the PDF. 
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From these examples, we learn that the use of the modified showering algorithm does not 
change the jet properties in a major way. Figure 7 compares the differential jet transverse 
energy distribution in Run I to CDF data. Here, jets are defined with R = 0.4 and are 
smeared with an energy resolution function AEt/Et = X.lj \JEx (in GeV). We consider 
the agreement between theory and data to be further evidence that our NLO corrections 
reproduce jet properties. 

VII. CONCLUSIONS 

We have presented a modified, parton showering algorithm that produces the total cross 
section and the event shapes beyond the leading order. These modifications are based 
on the Qr^space resummation. The parton showering itself is modified by using a new 
PDF (called mPDF) which encodes some information about the hard scattering process. 
Simultaneously, the explicit, hard emission is included, but only after subtracting out the 
contribution already generated by the showering: this correction is called Y . The presence of 
Y yields a smooth transition from the parton showering to single, hard emission. We modified 
the PYTHIA Monte Carlo to account for these corrections, and presented comparisons with 
Run I W boson data and predictions for W and Higgs boson production at the Tevatron in 
Run II. 

The scheme works very well at NLO for the cases considered in this study, and the 
correct cross sections, transverse momentum distributions, and jet properties are generated. 
We have compared our kinematic distributions to the case when the results of the standard 
showering are multiplied by a constant i^-factor to reproduce the NLO cross section. We 
find variations on the order of 10% for small and medium transverse momentum. 

There are several effects which still need study. We have not included the exact dis- 



tributions for the decay of the leptons |3^ for W and Z production, which are resummed 
differently. It is straightforward to include such effects. In the theoretical discussion and 
numerical results, we have focussed on initial state radiation, but our results should apply 
equally well for final state radiation. The situation is certainly simpler, since final state 
radiation does not require detailed knowledge of the fragmentation functions. Also, the case 
when color fiows from the initial state to the final state requires study. A resummed calcula- 
tion already exists for the case of deep inelastic scattering ||3^, and much theoretical progress 



has been made for heavy quark production ||38|. We believe that the modified showering 



scheme outlined in this study generalizes beyond NLO, just as the analytic calculations can 
be calculated to any given order. For example, we could include hard PF + 2 jet corrections 
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to Y. For consistency, however, higher order terms {A and B) may also need to be 
included in the Sudakov form factor. 

The modified PYTHIA subroutines used in this study and an explanation of how to use 
them are available at the web address moose.ucdavis.edu/mrenna/shower.html. 
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FIG. 1. The prediction of the W boson transverse momentum distribution in Run I at the 
Tevatron (sohd hne) compared to the D0 data (the upper and lower dashed hues represent the 
errors on the middle line). The prediction includes the effects of the modified parton distribution 
functions, the correction to the hard scattering process, and a tuned primordial k^ of 3.6 GeV. 
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W Boson Production at the Tevatron, Run II 
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FIG. 2. Predicted W boson transverse momentum distribution in Run II (solid line) showing the 
individual contributions from showering (long dashes), the corrected hard emission (short dashes) 
and the standard PYTHIA prediction using CTEQ4L structure functions (dot-dash). 
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100 GeV Higgs Boson Production at the Tevatron, Run II 
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FIG. 3. Same as Fig. 2, but for Higgs boson production {rriH = 100 GeV). 
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100 GeV Higgs Boson Production at the Tevatron, Run II 
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FIG. 4. The ratio of Higgs boson transverse momentum distributions {niH = 100 GeV) from 
the modified showering algorithm and from the standard showering algorithm using the CTEQ4L 
PDF. The CTEQ4L result has been multiplied by a constant i^-factor to reproduce the NLO rate. 
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FIG. 5. Jet multiplicity in W boson events for different Qt binnings. 
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FIG. 6. Jet multiplicity in Higgs boson events for different Qt binnings. 
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FIG. 7. Comparison of the differential jet cross section in W boson + jet events in Run I (dashed 
hne) to the CDF data (circles). The jets have been smeared by the resolution AEt/Et = 1.2/ ^/Et 
(in GeV). 
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